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To assist in the search for phenomenologically realistic models in the string landscape, we must develop tools 
for investigating all gauge charges, including {7(1) charges, in string models. We introduce the process for 
constructing fermionic string models and present an algorithm for determining the U(l) gauge states and 
U(l) charges in weakly-coupled free fermionic heterotic string (WCFFHS) models given their matter and 
non-Abelian gauge content. We determine the computational complexity of this algorithm and discuss a 
particular implementation that can be used in conjunction with a framework developed at Baylor University 
for building WCFFHS models. We also present preliminary results regarding matter state uniqueness for a 
run of 1.4 million gauge models and find that U(l) charges affect the number of unique matter states in a 
significant minority of models constructed. We conclude by considering future avenues of investigation to 
which this algorithm may be applied. 
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I. INTRODUCTION 

String theory models have the potential to completely 
describe both the particle content and field interactions 
of our universe, yet attempts to obtain precise, testable 
predictions based on string theory have been stymied by 
the sheer number of viable models. Work by a number 
of theorists in the early 2000 's helped establish that the 
number of stable models, though finite, could be upwards 
of 10 100 fir— ]. Any predictions based on string theory will 
require a more complete understanding of this collection 
of models, known as the string theory landscape. In par- 
ticular, systematic searches of the landscape can help us 
gain a statistical understanding of the range of available 
models. Such searches, however, require fast algorithms 
for determining properties of string models. In this pa- 
per, we discuss a new algorithm used to determine the 
U(l) charges of particles in weakly coupled free fermionic 
heterotic string (WCFFHS) models. Such models are 
of particular interest due to a number of studies which 
have discovered phenomenologically realistic models of 
this typei^. 



A. WCFFHS model building 

For our investigation of U(l) charges, we built on re- 
search conducted at Baylor University to develop a fast, 



versatile framework for building WCFFHS models 39 . 
The process of model-building begins with a set of 64- 
component basis vectors which specify information about 
the boundary conditions for fermions on the world-sheet 
as well as a valid GSO (Gliozzi, Scherk, and Olive) pro- 
jection matrix, which is used to eliminate non-physical 
particle states encountered during model-building. The 
64 components of the basis vectors are divided into two 
parts, the first of which specifies information about the 
left-movers on the world sheet and the second of which 
specifics information about the right-movers. The num- 
ber of components in each part is determined by the 
number of large spacetime dimensions. For four large 
spacetime dimensions, the left-moving part consists of 
20 components while the right-moving part consists of 
44 components. The 64 components of the entire vector 
are usually labeled as shown: 



,16\ 
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Additionally, each input vector has an order N, which 
specifies the denominator of all components in the basis 
vector. These components take on values between — 1 
and 1. 

These vectors must obey a number of constraints to 
constitute valid basis vectors. One particularly useful 
constraint specifies that every component in the entire 
vector must be paired with an equal component in order 
for the vector to be valid. This allows us to represent 
these vectors in a complex basis using only half as many 
components, an important advantage for manipulating 
them in calculations. In the description of our algorithm, 
we will use a complex basis to describe all vectors. 

The input basis vectors must also obey certain con- 
straints to guarantee that the generated model satisfies 
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the condition of modular invariance. This condition en- 
sures that modular transformations on the world-sheet 
will not alter the physics of the model. The specific con- 
straints on the input basis vectors {af } with orders {Ni} 
are given by 

Ni (aff = (mod 16) for even N (2) 

N% (af) 2 = (mod 8) for odd N (3) 

Nijdf ■ af = (mod 8), (4) 

where iVy is the least common multiple of Ni and Nj. 

Given a set of valid basis vectors, we construct sectors 
{di} by taking all linear combinations of the form 

a = rnjaf (5) 

where the rrij obey 

mj GZn[0,JVj). (6) 

The sectors determine how the free world-sheet fermions 
fi transform when transported around a non-contractible 
loop as given by 

/ — ► e -** a f. (7) 

This information is determined by the topology of the 
world-sheet. All possible charges or particle states Q a 
may then be constructed from these sectors as given by 

Q S = \a + F, (8) 

where the components of the fermion number operator F 
take on values given by 

F p G {1,0,-1}. (9) 

We maintain only states that are massless at the string 
scale, because these are the states that exist at low en- 
ergies and therefore represent observable particles. We 
exclude string-scale massive particles because they will 
not be observable. Additionally, we apply the GSO pro- 
jection to remove non-physical states such as tachyons. 
This projection ensures the remaining gauge states fit 
together to form a group theoretic root structure, com- 
patible with the remaining matter states 4 ^. The GSO 
coefficient matrix kij is subject to modular invariance 
constraints of its own, given by22 

Njkij=0 (mod 2) (10) 

kij + kji = -af ■ af (mod 2) (11) 

k u + ka = -af - Si (mod 2). (12) 

Given a valid GSO coefficient matrix, we reject any par- 
ticle states failing the requirement 

af -Q a = rnjhj + (mod 2) (13) 



where Sj equals if af is a bosonic state, and 1 if it is 
fermionic, and rrij and a are defined as in |Equation 5| 

We can identify the type of particle that a particular 
state specifies by considering the components of the left- 
moving part for that state. Specifically, for fermions 

^=^ 2 =±\, (14) 
and for gauge bosons 

V/ = ?A 2 = ±1, (15) 

with sign indicating chirality. Furthermore, all other left- 
moving components of gauge states are equal to zero. 
Thus, if the right-moving part of the U(l) gauge states 
can be generated, adding the left-moving part is straight- 
forward. Once all non-physical states have been removed 
from the model, we can determine the U(l) charges of 
the remaining physical matter states. 

B. 17(1) charges in WCFFHS models 

In order to fully understand the gauge interactions in 
a particular string model, we must determine the U(l) 
charges of the matter states. This point is perhaps 
best illustrated by considering the role of U(l) charges 
in the familiar context of the Standard Model. In the 
Standard Model, the charge of a particle under U(l) 
determines that particle's hypercharge, which in turn 
determines its electric charge through the Gell-Mann- 
Nishijima formula 41 

Q = h + j, (16) 

where Q is the electric charge of the particle, ^3 is an 
eigenvalue of its isospin, and Y is its hypercharge. 

In addition to completing the picture of the gauge in- 
teractions for a model, analyzing U(l) charges can help 
eliminate unwanted exotic particles from a model - in par- 
ticular fractionally charged exotic particles. If a certain 
type of U(l) charge known as an anomalous U(l) charge 
is present in a model, we can break that U(l), provid- 
ing a mass to the exotic particles at the string scale and 
thereby eliminating them from the low-energy effective 
field theory^ 5 . This can help dramatically in searching 
for phenomenologically realistic models, but identifying 
models of this type requires an efficient method for cal- 
culating U(l) charges. 

II. ALGORITHM FOR DETERMINING 17(1) CHARGES 

In order to determine the U(l) charges for the matter 
states in a string model, we must first construct the U(l) 
gauge states. These gauge states, by definition, have a 
number of convenient mathematical properties that sim- 
plify their construction. Firstly, because they are gauge 
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states, their left moving parts are already determined, as 
noted in Iscction J~K\ Secondly, we know that the right- 
moving part of the U(i) states must be orthogonal to the 
right-moving parts of the simple roots of the non-Abelian 
gauge groups in the model. If the state were not orthog- 
onal to the simple roots of any one of the non-Abelian 
gauge groups, it would be a part of that gauge group 
and not the U(l) gauge group. Finally, since represen- 
tations of U(l) have only a single generator, every U(l) 
gauge group in the model has exactly one gauge state and 
hence, the right-moving parts of every £7(1) gauge state 
must be orthogonal to the right-moving part of every 
other £7(1) gauge state. Thus, if we can construct a set 
of vectors which are orthogonal to the right- moving parts 
of the simple roots of the non-Abelian gauge groups in a 
model as well as to each other, we will have constructed 
the right-moving parts of all £7(1) gauge states in the 
model. Once the gauge states have been constructed, we 
take the dot product of the right moving part of each of 
the matter states with the right-moving part of the £7(1) 
gauge states to determine its charge under each £7(1 ) in 
the model. 

Throughout our discussion of the algorithm for 
generating £7(1) gauge states, we will use a complex 
basis for all vectors, as described in Isection J~K\ So, we 
can represent every pair of equal components in the 
real basis, using a single component in the complex 
basis. To simplify our discussion of the problem, we 
denote the set of vectors given by the right-moving parts 
of the simple roots of the non-Abelian gauge groups 
by ~Vsr- We denote the vectors whose components 
are given by the right-moving parts of the £7(1) gauge 
states by Vj/m- Furthermore, we denote the cardinal- 
ity of Vsr by Nsr and the cardinality of ~Vu(i) by Nym . 



be done by constructing a matrix whose rows correspond 
to the components of the vectors in Vsr and performing 
Gauss- Jordan elimination on that matrix. After elimi- 
nation, we identify all columns containing only O's, and 
we generate vectors in Vr/m which consist of all O's ex- 
cept for the component corresponding to one of these 
columns. Gauge states of this form are known as ex- 
ternal gauge states. We then proceed to generate any 
remaining vectors of Vj/m. 

Let nu(i) denote the number of vectors in Vym which 
have already been generated, where njjru may range 



from to N, 



Co- 



in order to construct an additional 



vector v of Vj/m from Vsr and the nj/m previously 
generated vectors, we must ensure that the dot prod- 
uct of v with every vector in Vsr, as well as with ev- 
ery previously generated vector in Vj/(x) is equal to 0. 
Since there are Drm components, this gives us a system 
of N$r, + nxjru linear constraint equations on Drm un- 
knowns. This means that the system has / free variables, 
where / is given by 



/ = D RM - Nsr - nu(i). 



(19) 



Thus we, can assign arbitrary values to / components 
of v and then straightforwardly solve for the remaining 
components of v by solving this system of equations. If 
we let F denote the value of / when 7i[/(i) equals 0, this 
further implies that the vectors of ~Vu(i) wm be unique 
only up to scaling and rotations in F dimensions. 

Although we must assign arbitrary values to / com- 
ponents of v, we are not entirely free in choosing the 
components to which we will be assigning values. This is 
perhaps most easily illustrated by considering the prob- 
lem of finding a vector w orthogonal to the following 
(3-dimensional) vectors: 



A. Outline of the algorithm 

Our algorithm for generating £7(1) gauge states begins 
with the set ~Vsr for a particular model. The dimension 
Drm of vectors in this set is given by 



D RM = 26 — -D, 



(17) 



where D denotes the number of large spacetime dimen- 
sions for the models. Since the vectors of V sr are 
formed from the simple roots of gauge groups, they are 
by definition linearly independent. Furthermore, since 
V(j(x) contains all vectors orthogonal to the vectors of 
Vsr and to each other, Vsr U Vjjm forms a complete 
basis. Therefore, 



N, 



(7(1) 



D 



RM 



Nsr- 



(18) 



If Drm — Nsr, then the model in question contains no 
£7(1) gauge groups. 

We start by attempting to generate any vectors in 
V;/(i) whose components all have value except for one 
component which takes on the value 1. This can easily 



1 2 2 
3' 3' 3 

2 1 1 
3' 3'3 



(20) 



Since there are two vectors given, the system of con- 
straint equations will consist of two equations in three 
unknowns. Thus, we may arbitrarily assign a value to 
one component of w. Let Wi denote the ith component 
of w, where i ranges from 1 to 3. If we choose to assign a 
value of 1 to w\, then (taking the dot product of w with 
the given vectors), the constraint equations are 



\w 2 + §«* 
3 3 



(21) 



Attempting to solve for w 2 and 103, we see that the con- 
straint equations reduce to: 



= 3 



(22) 
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The system is clearly inconsistent. A quick inspection of 
the original vectors shows that this occurred because they 
were not linearly independent in the subspace formed by 
the components for which we were solving. 

In order to avoid this problem in generating a vector 
v in Vj/(i) , we construct a matrix whose rows are deter- 
mined by the components of the vectors in *Vsr &s well 
as the vectors in ~Vu(i) which have already been gener- 
ated. Then, we perform Gaussian elimination on that 
matrix. After Gaussian elimination, we identify those 
columns which do not contain pivots and assign values 
to the corresponding components of v. We can safely do 
this since the vectors are now guaranteed to form linearly 
independent subspaces in the remaining components. 

To illustrate this method, consider the previous exam- 
ple of generating the 3-dimensional vector w. This time, 
we begin by creating a matrix from the given vectors: 



12 2 

I ! I 

3 3 3 



After performing Gaussian elimination, this gives 



2 1 1 

3 3 



(23) 



(24) 



Since the third column of this matrix does not contain a 
pivot, we can safely assign a value (arbitrarily chosen to 
be 1) to W3. Proceeding as before, the system of equa- 
tions for the remaining components is given by 

w\ + 2w2 = —2 

2w x +w 2 = -1. (25) 
This system is consistent (as expected), and we get: 



w = (0,-1,1) 



(26) 



After assigning all necessary values to the components 
of v and solving for the remaining components, we add 
v to the set of generated vectors and repeat until we 
have produced the A r c/ ( 1 ) vectors in ~Vu(i) ■ These vectors 
specify the components of the right-moving parts of the 
U(l) gauge states, and we can trivially prepend the 
necessary left-moving part to make the full state a gauge 
state. Then, taking the dot product of the right-moving 
part of each U(l) gauge state with the right-moving 
part of each matter state in the model, we produce the 
charges of each matter state under each U(l) gauge 
group in the model. 



B. Implementation details and additional considerations 

Our implementation of the algorithm was written as an 
extension of a framework developed at Baylor University 
to build WCFFHS models 39 . Given a set of valid input 
vectors as well as a valid GSO coefficient matrix, this 
framework generates the matter and non-Abelian gauge 



content for the model. Our algorithm then generates 
the U(l) charges for the generated matter states. The 
entire framework, including our implementation of the 
algorithm, is written in CH — h, and care is taken through- 
out to avoid floating numeric types in order to eliminate 
rounding errors. The framework was written with exten- 
sibility and readability in mind, since we hope that other 
researchers will eventually be able to use this software to 
further explore the string landscape. In order to guide 
future development, a few choices in the details of our 
implementation deserve particular consideration. 

In our implementation, Gauss- Jordan elimination is 
used to solve for the unassigned components in the gen- 
erated vectors. Since Gaussian or Gauss- Jordan elimi- 
nation is required for two other steps in the algorithm, 
this allows for significant code re-use, but a few down- 
sides to Gauss- Jordan elimination must also be consid- 
ered. Firstly, Gaussian elimination is not a numerically 
stable algorithm, so we must consider the possibility that 
the program would overflow the CH — h long integer type 
used to store the coefficients of the systems of equations. 
In order to minimize this possibility, we implemented 
partial pivoting, which substantially improves the sta- 
bility of Gaussian elimination's,. We then tested the 
program by asking it to produce a complete orthogonal 
basis from a single randomly generated 22-dimensional 
vector. This goes beyond the worst case that our imple- 
mentation would attempt to solve, since all non-Abelian 
gauge groups that could appear in WCFFHS models have 
at least two generators. After running several hundred 
thousand trials, we determined that so long as the abso- 
lute value of the components in the randomly generated 
vector remained below about 250 (a very reasonable as- 
sumption) , the program never overflowed the long integer 
type. Furthermore, since this test went well beyond the 
worst case that our implementation will reasonably en- 
counter, the absolute value of the components of the in- 
put vectors can likely range substantially higher without 
difficulty. 

The second downside to Gauss- Jordan elimination is 
that it is typically somewhat slower than other common 
methods for solving systems of linear equations' 2 . This 
downside is balanced by a number of factors, including 
the fact that Gauss- Jordan elimination with pivoting is 
usually slightly more stable than most other general so- 
lution methods, and the fact that Gaussian or Gauss- 
Jordan elimination is required in other steps of the algo- 
rithm. Furthermore, since major goals of our implemen- 
tation included readability and extensibility, the use of 
a single easily-understood algorithm was deemed prefer- 
able to several faster, more esoteric algorithms. Never- 
theless, if future advances improve the speed of generat- 
ing the matter and non-Abelian gauge content of string 
models to the point that the process of generating U(l) 
charges becomes a bottleneck, other solution methods 
should be considered. 
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III. ANALYSIS OF COMPUTATIONAL COMPLEXITY 

In order to facilitate analysis of the computational 
complexity of the algorithm for generating all U(l) gauge 
states, we consider each significant step in the process. 
Using the notation of Isection III this process begins with 
Gauss- Jordan elimination on an N$r x Drm matrix to 
find any external U(l) gauge states. Gauss- Jordan elim- 
ination is generally O (n 2 m) on the number of rows (n) 
and columns (m) in a matrix, but our implementation 
includes a sorting step after each elimination step in or- 
der to maintain optimal pivoting. At worst, this sort 
will be an O (DrmN$r 2 ) operation, but is likely to be 
closer to or better than the average case, which is an 
O (DjiMNsRlogNsR) operation. Thus, the entire oper- 
ation of Gauss- Jordan elimination with pivoting will have 
complexity 

O (D RM N SR 3 ) (27) 

for the worst case. 

Next, we must generate all internal U(l) gauge states. 
If we let Next denote the number of external 17(1) gauge 
states, then there will be Drm — N$r — N eK t such states. 
For each state, we must perform Gaussian elimination to 
determine to which components we may assign arbitrary 
values. We must also perform Gauss- Jordan elimination 
to determine the values of the remaining components. In 
the first of these steps, the sorting process is simplified by 
the fact that the sorted results from the previous iteration 
can be stored and used to analyze the next vector to 
be constructed. Thus, the worst case for the Gaussian 
elimination step will be: 

0(D RM (NsR+n u{1) ) 2 y (28) 

where nj/m is the number of previously generated U(l) 
gauge states, including the external U(l) gauge states. 
The Gauss- Jordan elimination step will have no such ad- 
vantage. Thus, using the same analysis as before, the 
worst case complexity for this step will be 

O {{Drm - (N SR + n u{1) )) (N SR + n£/(i)) 3 ) . (29) 

Both the Gaussian and Gauss- Jordan elimination steps 
will be repeated until Jiwi) = Drm — Nsr- The com- 
plexity of this step will be at a maximum when 

N S R + n u{1 ) = \DRM. (30) 

Thus the overall worst case complexity of the algorithm 
for generating {7(1) gauge states will be 

O (D RM 4 ) . (31) 

Therefore, the overall worst case complexity depends 
solely on the number of right-moving components in the 
states of the model, which is itself determined by the 
number of large spacetime dimensions. 
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FIG. 1. A linear regression of the number of Z7(l) groups vs 
the number of SU(4) groups for models in our initial run 



IV. INITIAL RESULTS 

In order to develop an initial understanding of the im- 
pact of U(l) charges on analysis of WCFFHS models, we 
constructed and analyzed about 1.4 million models and 
their U(l) charges. This was done using the framework 
developed at Baylor for building WCFFHS models and 
our implementation of the algorithm for generating U(l) 
charges. For this run, we used only non-supersymmetric 
layer-1 gauge models with an order of 22 or less. Fur- 
thermore, we restricted our run to models which contain 
at least one copy of the Standard Model gauge group, as 
established by previous runs on model-building software 
developed at Baylor. While this is a fairly specific subset 
of the WCFFHS landscape, this run can help provide an 
initial sense of how U(l) charges affect our understand- 
ing of previous searches of the landscape. We present a 
few of the more interesting results of this run here. 

The models studied had anywhere from three to eight 
U(l) gauge groups, with a median value of four. Since 
the WCFFHS model-building framework produced only 
the non-Abelian content of its models prior to our re- 
search, we were interested if the presence or absence of 
any particular non-Abelian groups was correlated with a 
higher number of U(l) groups. Simple regression anal- 
ysis showed no significant linear correlation between the 
number of any particular non-Abelian group in the model 
and the number of U(l) groups in the model. In fact the 
highest R 2 value for a linear regression between the num- 
ber of occurrences of a particular group and the number 
of U(l) groups in the model was only 0.49. This result 
was for SU(4) groups versus U{\) groups as shown in 
|Figurc T| As this figure indicates, even the correlation 
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FIG. 2. Distribution of B for all models built 



between SU(4) and U(l) groups is almost entirely in- 
significant. 

Next, we examined the role of U(l) charges in deter- 
mining the uniqueness of matter states within a model. 
While examining uniqueness between models would be 
of interest, the fact that U(l) gauge states are unique 
only up to rotation prevents direct comparison of U(l) 
charges, and related problems affect the non-Abelian 
charges- 9 . Thus, we confine our analysis to compar- 
isons of matter states within a model. In particular, we 
were interested in the ratio of matter states unique in 
their non-Abelian charges to matter states which were 
truly unique (i.e. unique in all charges, including U(l) 
charges). We refer to this ratio as B, and it gives us 
some sense of how great of an error would be made in 
studies which only examined non-Abelian charges when 
attempting to identify unique matter states. The distri- 
bution of B for our run of 1.4 million models is shown in 
|Figure 2| As this figure indicates, consideration of J7(l) 
charges does not affect the number of unique particles for 
the majority of models in our study, since B had a value 
of one for most of these models. A sizable minority had 
values less than one, however. Therefore, we conducted 
linear regression analyses of B vs the occurrences of each 
group in our models in order to see if there were any 
correlations that could be used to predict the impact of 
U(l) charges on matter state uniqueness prior to building 
the charges. These analyses showed no significant linear 
correlation between the value of B and the occurrences 
of any given group [including U{1) groups] in the model. 
The highest R 2 value among these regressions was only 
0.07. 

After examining the effect of U(l) charges on matter 
state uniqueness, we considered the overall distribution of 



FIG. 3. Distribution of unique matter states for all models 
built 

the number of unique matter states for the models in the 
run, as shown in |Figure~3| As the figure indicates, this 
distribution is roughly symmetrical with a sharp peak at 
12 unique matter states. The mean number of unique 
matter states was 12.27, with a standard deviation of 
2.43. A Shapiro- Wilks normality test on these data gives 
sufficient evidence to reject the hypothesis that they fol- 
low a Gaussian distribution {p < 0.01). It is possible, 
however that this lack of normality is a feature only of 
the specific population chosen for this study and does not 
hold more generally for the WCFFHS landscape. While 
all of these results are preliminary, they give us some 
sense of how the construction of U(l) charges can affect 
our understanding of matter state uniqueness for models 
in the landscape. 



V. CONCLUSIONS 

In this paper, we presented a new algorithm for cal- 
culating U(l) charges of particles in WCFFHS models 
given the matter and non-Abelian gauge content for those 
models. We analyzed the computational complexity, and 
found it to be O (Drm 4 ) in the worst case, where Drm 
refers to the number of right-moving components of the 
particle states in a complex basis. We also discussed 
an implementation of this algorithm developed as an ex- 
tension of a model-building framework developed at Bay- 
lor University which emphasizes readability, extensibility, 
and speech, and we presented results on matter state 
uniqueness from an initial run of 1.4 million gauge mod- 
els using this implementation. Specifically, we found that 
the U{1) content affected the number of unique matter 
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states only in a minority of models. We found no strong 
correlations that would help predict which models lie in 
this minority, however. 

Determining the U(l) charges for a particular model 
is important in understanding its gauge content, but it 
can also help in constructing phenomenologically realis- 
tic models by identifying anomalous U(l) charges. Addi- 
tionally, constructing U(l) charges may shed some light 
on the problem of identifying unique models, since mul- 
tiple inputs may all lead to models with the same matter 
and non-Abelian gauge content. While some studies have 
been conducted to analyze patterns in the occurrence of 
unique models, it is possible that by studying the U(l) 
gauge content, we will find that some models with the 
same matter and non-Abelian gauge content are actually 
distinct due to their varying U(l) charges^. 

Like many aspects of the string landscape, we cannot 
predict everything that the U(l) gauge content of string 
models will tell us. Only through systematic exploration 
of the landscape can we completely understand how U(l) 
charges affect the search for phenomenologically realistic 
models. 
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